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Abstract 

The finite-size effects in critical phenomena of a thin film system are stud- 
ied from a mean field (MF) approach with <f> 4 model for second-order phase 
transition. The influence of boundary condition on the critical properties are 
emphasized. Scaling functions for relative free-energy and relative specific 
heat capacity are given. 

PACS numbers: 68.35.Rh, 05.70.Jk 

The critical phenomena in the thermodynamic limit is characterized by the di- 
vergence of the correlation length £ near critical point. Nowadays, the experimental 
techniques have become so advanced that the correlation length £ can be pushed 
up to several thousand A and the samples under study become comparable with £. 
As a consequence, the effects of finite size of the samples on the critical phenomena 
become increasingly important. Generally such effects depend on the shape of the 
sample, on the boundary condition, on the dimension of the system and on the num- 
ber of components of the order parameter. Finite-size effects near critical point have 
been remained to be a topic of active research over the past two decades. In fully 



finite or quasi-one- dimensional systems the phase transition is smeared out, whereas 
in thin films of thickness L the critical temperature Tc(L) is shifted with respect 
to the bulk Tc- Without any doubt, detailed theoretical studies are necessary in 
order to understand the most important finite-size scaling structures predicted by 
phenomenological [1] and renormalization-group theories [2]. 

Since 1985 various perturbative approaches to calculations of finite-size scal- 
ing functions with the 4 model have been developed [3-6]. All those perturbative 
approaches are based on the division of the free energy into a Gaussian term and 
higher order perturbative terms. In fact, the coupling constant with the 4 model of 
the perturbation is not necessary small, so that the convergence of the perturbation 
expansion cannot be ensured. Thus some more effective approaches to the calcu- 
lation are needed. It is clear that the random changes of the order parameter due 
to thermal motion should be around certain average and the fluctuations around 
the average are generally not too large except in the temperature region close to 
the critical one. If the averaged order parameter can be obtained, the fluctuations 
can be described as a Gaussian term and small higher order terms. MF calcula- 
tions are important. First of all, they furnish the lowest order approximation of the 
real problem, and some characteristics can be obtained from such calculations. The 
advantage of the MF approximation is that it allows a systematic discussion of the 
corrections to that approximation. For an infinite size system, such a MF calculation 
is quite easy and can show some features of critical phenomena. If the same results 
are used for system with finite size, thermodynamic quantities will have the same 
jumps as the bulk system, thus are physically unacceptable for finite size systems. 
The normal solution to that problem is with the zero-mode approximation instead 
of the MF approach. The fundamental problem in the usual MF approach when 



used to finite size systems, it seems to us, is that there exist boundary influences. 
For system with boundaries, one should consider the influences of the boundaries on 
the thermal properties near the bulk critical point. In addition, the spatial distribu- 
tion of the order parameter should be taken into account for finite size system. In 
the MF calculations for an infinite system, the gradient term has no contribution, 
because the system prefers a uniform distribution of order parameter. For a finite 
size system, however, the spatial distribution cannot be considered as uniform any 
longer due to the influence of the boundary though the condition of minimum free 
energy would prefer a smooth distribution. Thus the gradient term should not be 
ignored in the MF calculations of finite size system. When this term is taken into 
account for a finite size system, the MF calculations become much more difficult 
than those for the infinite system. On the other hand, MF calculations provide 
the base for the consideration of fluctuations. The MF approximation is just the 
tree level approximation of the corresponding field theory, and corrections to mean 
field can be given by suitable perturbation theories. By the way, the calculations 
can be expected to give reasonable results and are very easy comparing to those in 
perturbation expansion. Exact calculations show that the MF approximation be- 
comes exact for an infinite system when the dimension d of the system becomes very 
large. To the best of our knowledge, no such work for finite size system has ever 
been done yet. In this paper, we try to show some features of finite-size effects in 
second-order phase transition from a MF approach suitable for limited systems. For 
simplicity, we only consider 4 model for a thin film under second-order phase tran- 
sition with a one-component order parameter as an example to show the method. 
The extension to fully finite system and to the case of n-component order parameter 
is straightforward. 



In Ginzburg-Landau theory of phase transition, the partition function for second- 
order phase transition is given as 

Z = Jv<j> exp(-if) (1) 

with effective Hamiltonian 



H = I d d x 



|0 2 + ^(V0) 2 + ^0 4 



(2) 



with 7 = a(T-T c ), a > 0. 

In the following, we consider a thin film with thickness L and assume the bound- 
ary condition for the system is of Dirichlet type for simplicity. For non-homogeneous 
boundary conditions, similar calculations can be done. The MF is the configuration 
of which makes the effective Hamiltonian H minimum, or simply 5H = 0. From 
the point view of partition function, MF theory is a saddle point approximation. 
The corrections to the approximation can be considered by expanding around the 
saddle point. From the condition of minimum of H, one gets 

V 2 = 70+^0 3 (3) 
with boundary condition 0(0) = 0(L) = 0. Since one intends to choose the config- 
uration of with H minimum, the averaged order parameter should vary only in 
the direction of L which will be chosen as x direction in later discussions. If 7 > 0, 
last equation has only trivial solution = 0, the same as for the bulk system. So one 



needs only to discuss the case with 7 < 0. Defining O = ^—67/6 and £ 2 = — 1/7, 
the last equation can be normalized as 

= -y + *3 , (4 ) 

with ^ = 0/0o and \&(0) = ^(L) = 0. Obviously, £ is the usual correlation length 
and is MF relative to 0o, the MF for the corresponding infinite system. Appar- 
ently, Eq. (4) is invariant under ^ <-> — ^ since ±^f correspond to the same effective 



Hamiltonian H. Then one can assume \I/ > in later discussions. Last equation 
can be integrated as 



with c an integral constant which, as will be shown afterward, is in the range [0,1/2]. 
It should be kept in mind that the solution of last equation should make the effective 
Hamiltonian H minimum. Thus the spatial change of ^ should be as slow as possible. 
Due to the symmetric boundary condition, \& will reach its unique maximum at 
x = L/2, where d^/dx = 0. Let k 2 = (1 — c — \/l — 2c) /c, the integral constant c 
as a function of £ and L can be determined through 



where F(k) is the first kind of complete elliptic integrals. From above expression, 
one can see that the integral constant c depends not on £ and L separately but on 
the combination of L/£. Thus it is convenient to define a reduced length of the 



system I = L/£ = Z y 1 — T/T C - The value of l depends only on the nature of 
the substance under study. In fact both I and l are dimension-less, they are ratios 
between the real thickness of the film and certain characteristic lengths. In fact, 
Eq. [4] can be written in terms of dimensionless variable I. Since k should be in the 
interval [0, 1] the integral constant c should be in [0, 1/2]. Smaller c corresponds 
to smaller reduced length /. The behavior of c with / is shown in Fig. I. A special 
feature appears as shown in the figure. If the reduced length I < ir, c = thus the 
maximum of the MF or averaged order parameter ^ = 0. In other words, under such 
condition the system will remain in the disordered phase despite the bulk critical 
point has been crossed from above. One should notice that what plays a role is not 
the real length L of the system but the reduced one / which depends also on the 




(5) 



L = 2£VTTk 2 ~F{k) , 



(6) 




deviation from the bulk critical point. This shows that if the reduced length of the 
system is too short, either due the small size of the system or due to its closeness 
to the bulk critical point, the influence of the boundary will dominate. This result 
is not surprising considering the physical meaning of the correlation length £ which 
can be regarded as an indication of the influence range of the boundary condition. 
The minimum values of the reduced length l min = ir, as a mathematical result, 
can be understood physically. If c is very small, the maximum of \I/ is also very 
small, so the \l/ 4 term in Eq. (||) can be neglected. Then the solution of Eq. (|J) is 
^ = y / csin(a;/£). As discussed above, \l/ should reach its maximum before x = L/2, 
thus I = L/£ > 7r, otherwise \1/ = 0. Of course, the existence of \1/ 4 term will modify 
the behavior of \l/ and make the two parts for x > L/2 and x < L/2 join together 
smoothly for larger I. 

The existence of a minimum reduced length l m [ n is an indication of the shift of 
critical point for a the finite size system with respect to the bulk one. One can 
easily get the critical temperature Tc(/ ) at which the finite system begins to have 
non-zero averaged order parameter 



If the characteristic length 1$ of the system goes to infinity (thus the system becomes 
unlimited), Tq{Io) in last equation returns to the bulk value Tq. But for small Iq the 
shift of critical temperature will be important. Specially, if Iq — > it, Tc{Iq) — > 0, 
so that the limited system will stay in the disordered phase until the temperature 
becomes absolute zero. Thus there exists a lower bound for the size of the system 
above which there may exist phase transition. With the increase of Z, the constant 
c increases monotonously and saturates at 1/2 for I — > oo. When c ~ 1/2, the MF 




(7) 



cf) almost equals to the value 0o except in the narrow bands near the boundaries. 
This shows that the effects of the boundary condition can be ignored for system 
with the reduced length I » 1. This point can be seen from the behaviors of other 
thermodynamic quantities. 

The effects of the boundaries on the distribution of MF \1/ can also be shown 
from Eq. (|5]). One can first normalize the film by variable transformation x = x'L. 
Then 

x> = iv / T+PF(sin- 1 ^ V ^\ k) , (8) 

where F(8, k) is the first kind of incomplete elliptic integral. From last expression 
the MF \l/ can be solved as a function of x' and Due to the existence of minimum 
reduced length Z min = 7r one is interested in the varying behavior of the distribu- 
tion of ^ verse V = I - ir. For /'=0.0625, 0.214, 0.328, 0.409, 0.863 and 10.00 the 
distributions of MF \1/ = ^(x',l) are shown in Fig. 2. From those figures one can 
see that the maximum value of \1/ increases with the increase of I and when I is 
large enough there is a plateau in the spatial distribution of The larger /, the 
larger the fraction of the plateau, indicating the lesser influence of the boundary 
condition. Because of the continuous change of the distribution with the reduced 
length, the singularities accompanied with the infinite system are smeared, as de- 
manded from general physical consideration. This can be shown from the behaviors 
of thermodynamical quantities. 

To study the finite-size effects on free energy, define = — 37 2 /26, which is 
the density of free-energy for infinite system. Then the ratio F r between free energy 
of the thin film per unit thickness F/L and that of the infinite counterpart per unit 



volume Foo can be shown as 



F 2 r l 

F r = — — = - / c(t)dt . (9) 



LFoo I Jo 

This ratio depends only on the reduced length of the system, indicating the mutual 
influences of the boundaries and intrinsic thermal motions. From the saturation of 
c(l), it is obvious that F T also saturates to 1 for I — > oo. But as shown in Fig. 3, the 
saturation of F r to its maximum value 1 is very slow, so to speak, the influences of 
the boundary condition is in force even the length scale L of the system is over ten 
times larger than the correlation length £. Thus the influences of boundaries should 
be considered even for systems with size scale much larger than £. 

Once the free-energy of the limited system is known, other thermal quantities 
can be obtained. As an example, one can consider the specific heat capacity. Define 
C r the ratio between specific heat capacity of the limited system and that of the 
infinite counterpart. After simple algebras, one arrives at 

c-OvWs + | Fr + L£. (10) 

CV(oo) 4 8 4 dl 
Again this ratio C r only depends on the reduced length of the system. The vary 
of C r with the reduced length / are shown in Fig. 4. In the region of / close to its 
critical value n, the relative specific heat capacity increases quickly with /, but the 
rate of increase slows down and finally saturates to 1 as expected. 

As a summary, we have investigated the finite-size effects on critical behaviors 
based on a MF approach adapted for the limited system with a one-component 
order parameter. Though the fluctuations have not been considered, interesting 
critical phenomena due to finite size of the system have been shown. The relative 
thermodynamical quantities such as the free energy and specific heat capacity with 
respect to those for infinite counterpart are expressed as scaling functions of the 



reduced length I. It will be fruitful to study the modification to the MF results from 
fluctuations, which will be discussed elsewhere. 
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Figure Captions 

Fig. 1 The behavior of the integral constant c as a function of the reduced length 
/ of the system. 

Fig. 2 Distributions of mean field \l/ as defined in text over the normalized space 
x' for various I' = I — n. 

Fig. 3 The ratio between the free energy density of a limited system and that of 
the infinite counterpart as a function of the reduced length of the system. 



. 4 The ratio between the specific heat capacity of a limited system and that 
of the infinite counterpart as a function of the reduced length of the system. 



